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Abstract: We generalize the local model of primordial non-Gaussianity by promoting the 
parameter /nl to a general scale-dependent function /nl(^)- We calculate the resulting 
bispectrum and the effect on the bias of dark matter halos, and thus the extent to which 
fNh(k) can be measured from the large-scale structure observations. By calculating the 
principal components of /nl(^)j we identify scales where this form of non-Gaussianity is 
best constrained and estimate the overlap with previously studied local and equilateral 
non-Gaussian models. 
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1. Introduction 

Primordial non-Gaussianity provides cosmology one of the precious few connections be- 
tween primordial physics and the present-day universe. Standard inflationary theory, with 
a single slowly rolling scalar field, predicts that the spatial distribution of structures in 
the universe today is very nearly Gaussian random (e.g. [1, 2, 3, 4, 5]; for excellent recent 
reviews, see [6, 7]). Departures from Gaussianity, barring contamination from systematic 
errors or late-time non-Gaussianity due to secondary processes, would be a violation of this 
standard inflationary assumption. Constraining or detecting primordial non-Gaussianity 
is therefore an important basic test of the standard cosmological model. 

Most of the study of non-Gaussianity in the literature to date has been carried out 
assuming the magnitude of departure from Gaussianity is scale-independent (e.g. [8, 9, 10]). 
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However, the assumption that /nl is constant for a wide range of scales could be an over- 
simplification, since the primordial cosmic perturbations were presumably produced from 
the time-dependent dynamics in the early universe. In particular, single-field inflationary 
models with interactions, along with most multi-field models, generically produce scale- 
dependent non-Gaussianity. It is therefore not surprising that scale-dependence of non- 
Gaussianity has been discussed in the community in recent years [11, 12, 13, 14, 15, 16, 
17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29]. Notably, the parameterization of the 
scale-dependent non-Gaussianity in our analysis is applicable to the curvaton [30, 31, 32, 
33, 34] and the modulated reheating scenarios [35, 36], which are of great interest for their 
potentially observable scale-dependent non-Gaussianity 1 . 

Motivated by such inflationary models that predict detectable scale-dependent non- 
Gaussianity, as well as a desire to have an easily usable basis for studying those models, 
we present a novel scale-dependent ansatz for primordial non-Gaussianity: we promote the 
parameter /nl to a free function of wavenumber /nl(&)- We define our model (Sec. 2), 
predict clustering bias of dark matter halos in our model (Sec. 3), obtain an upper bound 
on the accuracy with which these new parameters could be measured with a future large- 
scale structure survey (Sec. 4), and compare our model with other parameterizations of 
non-Gaussianity in the literature (Sec. 5). 

2. Scale dependent non-Gaussianity 

The most commonly discussed model of non-Gaussianity, often referred to as the local 
model, is defined via [8] 



Here, $ denotes the primordial curvature perturbations (Bardeen's gauge-invariant po- 
tential), 4>g(x) is a Gaussian random field, and the constant /nl is the non-Gaussianity 
parameter. The local model has been much studied, in part because it is the first two 
terms of the most general local form of non-Gaussianity [40]. 
In Fourier space, Eq. (2.1) becomes 



(Hereafter, we omit the subscript G on the Gaussian distribution when it is clear from 
context.) In this paper, we study a model that generalizes Eq. (2.2) - we allow /nl to 
vary with k as well, while assuming isotropy and homogeneity (so /nl(&) = /nl(&))- The 

1 For instance, when the observed perturbations originate from the single curvaton field, the "running" 
(with scale) of the non-Gaussianity parameter is proportional to the third derivative of the curvaton po- 
tential, V'" [37, 38, 39]. Given that this third derivative is not tightly constrained from the observed 
power spectrum, it can potentially lead to observable and scale-dependent non-Gaussianity. Therefore, 
constraints on the running of non-Gaussianity can be a powerful probe of the origin of the primordial 
curvature perturbations. 
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gravitational potential in the new model is defined via 

/rfik' 
—ftk'Mk-k'). (2.3) 

As mentioned above, this form of non-Gaussianity is expected in curvaton or modulated 
reheating scenarios (see e.g. Ref. [37], where this form explicitly appears in the study of 
these models). 

Note that this new ansatz is not local, which is clear when we transform back into real 
space: 

*{x) = <t> + f NL (x) * (<^) 2 - <<Xx) 2 », (2.4) 

where * represents convolution and x denotes a three-dimensional spatial coordinate. These 
primordial perturbations $>(k) are related to the present-time (z=0) smoothed linear over- 
density 6r by the Poisson equation: 

2 k 2 T(k) 

<**(*) = o^2?rW R (k)m = M R (k)*(k); (2.5) 

O II q S l m 

where T(k) is the matter transfer function, Ho is the Hubble constant, £l m is the matter 
density relative to critical today, and Wn(k) is the Fourier transform of the top-hat filter 
with radius R. The smoothing spatial scale R is related to the smoothing mass scale M 
via 

M = ^R 3 Pmfi , (2.6) 

where p m <d is the matter energy density today. The choice of mass scale is discussed further 
in section 4.1. 

The bispectrum in our generalized model becomes 

B^h, k 2 , fc 3 ) = 2[/ NL (A ; i)P^(A; 2 )P^(fc3) + perm.], (2.7) 

where is the power spectrum of potential fluctuations. This reduces to the familiar 
expression B(k±, k 2 , k$) = 2 /nl (■?*</> (fci)-P</>(fc2) +perm.) when /nl is a constant. 

Notice the difference between our ansatz for the scale-dependent /nl(&) (which has 
the corresponding bispectrum Eq. (2.7)) and the particular form of scale-dependent non- 
Gaussianity, discussed elsewhere in the literature, which is defined as /nl(&i> ^2, ^3) = 
B^ki, k 2 , k^) /[2P ( p(ki)P ( f l (k2) + perm.] ([25, 27, 26]). The two forms are inequivalent, and 
either form can be borne out in realistic inflationary models; however, given that our form 
lives in a lower-dimensional A:-space, it is easier to simulate it numerically [41] or treat it 
with the Fisher matrix analysis, as we do in this paper. 



3. Non-Gaussianity and Bias 

3.1 The effect of a non- vanishing bispectrum on bias 

Dalai et al. [42] found, analytically and numerically, that the bias of dark matter halos 
acquires strong scale dependence if /nl 7^ 0: 

on tf2 

H*)=^ + /NL(6.-i)4 ag(a) ; w v (3.1) 
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Here, bo is the usual Gaussian bias (on large scales, where it is constant), 5 C ~ 1.686 
is the collapse threshold, a is the scale factor, $7 m is the matter density relative to the 
critical density, Hq is the Hubble constant, k is the wavenumber, T(k) is the transfer 
function, and g(a) is the growth suppression factor 2 . This result has been confirmed by 
other researchers using a variety of methods, including the peak-background split [43, 44, 
45, 46], perturbation theory [47, 48, 49], and numerical (N-body) simulations [50, 51, 52]. 
Astrophysical measurements of the scale dependence of the large-scale bias, using galaxy 
and quasar clustering as well as the cross-correlation between the galaxy density and CMB 
anisotropy, have recently been used to impose constraints on /nl already comparable to 
those from the cosmic microwave background (CMB) anisotropy [45, 43], giving /nl = 
28 ± 23 (lcr), with some dependence on the assumptions made in the analysis [45]. In the 
future, constraints on /nl are expected to be on the order of a few [42, 53, 54, 55]. The 
sensitivity of the large-scale bias to other models of primordial non-Gaussianity has not 
yet been investigated much (though see analyses in e.g. [56, 57]). 

Following the MLB formula [58, 59], one can express the two point correlation func- 
tion of dark matter halos, ^(xi, x 2 ), in terms of certain configurations of the correlation 
functions of the underlying density field, £ R N ^ ■ In the high-threshold limit [y 3> 1), this 
becomes: 



t,h(xi,x 2 ) 




00 iu 
V 



Xl,...,Xl, x 2 ,...,x 2 
j times (N — j) times 



(3.2) 



where Xij = \x,i — Xj\,u = 5 c /a R represents the peak height, and £R^ n \r) is the n-point 

correlation function of the underlying matter density smoothed with a top-hat filter of 

radius R. Keeping the terms up to the three-point correlation function, which would be 

reasonable for the observationally allowed range of /nl ; the expansion series gives us the 

halo correlation function in terms of the field correlation functions: 

2 3 
£h(x 12 ) = - T ^(x 1 ,x 2 ) + - T ^\x 1 ,x 1 ,x 2 ). (3.3) 

°R °R 

The Fourier transform of the real-space correlation function - the power spectrum - 
is given, to the same expansion order as Eq. (3.3), by 

P h (k) = ^P R (k) + ^ [ ^B R (k,q,\k-q\) + ... (3.4) 

a R a R J \ Z7T ) 

The first term on the right-hand side includes the familiar (Gaussian) bias b = vja R (in 
the high-peak limit for which the MLB formula is valid) for the Gaussian fluctuations. The 
effects of non-Gaussianity on the galaxy bias are represented by the second term, including 
the bispectrum Br, which vanishes for the Gaussian fluctuations. 



2 The usual linear growth D(a), normalized to be equal to a in the matter-dominated epoch, is related 
to the suppression factor g(a) via D(a) — ag(a), where g(a) is normalized to be equal to unity deep in the 
matter-dominated epoch. 
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3.2 Prom the bispectrum to bias 

If we denote the full bias of dark matter halos by b + Ab, where b represents the bias for 
the Gaussian fluctuations and Ab is the non-Gaussian correction, then 

2 



El 
PR 



Ab\ 

1 + " 



(3.5) 



where -f\ and Pr are the power spectra of halos and dark matter, respectively. The non- 
Gaussian correction to the linear peak bias to the leading order becomes 

d 3 q 



Ab 



(k) 



1 



or 2P R {k) J (2tt) 



rB R (k,q, \k - q\ 



(3.6) 



where Br is the matter bispectrum on scale R. Hence, the non-Gaussian correction Ab(k) 
can be expressed in terms of the primordial potential fluctuations as ([44]): 



Ab 



(k) 



S c 



D(z) 87r 2 a R M R (k) J 



dkik\MR{ki) 



l P^k) 



(3.7) 



We perform the integration over all triangles. The triangles' sides are k±, k 2 , and k; the 
cosine of the angle opposite k 2 is fi, so k 2 = kf + k 2 + 2k\k\x. Ai R {k) is the same function 
denned in Eq. (2.5), and the time dependence of the critical threshold for collapse is given 
as 5 c (z) = 5 C /D(z), with 5 C = 1.686. 

3.2.1 Constant /nl 

Eq. (3.7) leads to the famous scale-dependent bias formula in the case of a constant /nl- 
For this model, the bispectrum is 



B^kt, k 2 , k 3 ) = 2/ NL [P ( j>{k 1 )P^(k 2 ) + perm. 
Through Eq. (3.7), this leads to the result 



(3.8) 



Ab 



(k) 



2/i 



NL 



D(z) 87T 2 a 2 R M R (k) 

2/nl^c -F(fc) 
D(z) M R (kY 



dhkfMRik^P^h) / dfiM R (k 



p^M 

pm 



+ 2 



where 



2^-2 



iir^o 



dkikfMRik^P^ki) / d^M R {k 2 



P*{k 2 ) 
P<t>(k) 



+ 2 



(3.9) 



(3.10) 



Note that there is a factor of 2 in Eq. (3.9) because we can exchange the order of integration 
of terms corresponding to k\ and k 2 . 

Finally, we rewrite Eq. (3.9) by denning 

1 



?x{k) 



Mk) 



8ir 2 o 2 R M R (k)P^k) 



dhkiMRik^P^h) / dnM R (k 2 )P<i,(k 2 ) (3.11) 



87r 2 a 2 R M R (k) 



dhkfMRiki^ih) / dfiM R (k 2 ). 



(3.12) 
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Then, for constant /nl, 



Ab 



(*) 



b v y D(z) 
and the derivative with respect to /nl is 



[^i(fc) + ^2(*;)], 







Ab 



(k) 



25 c 
D(z) 



(3.13) 



(3.14) 



<9/nl 

3.2.2 Scale-dependent /nl 

Now we repeat the analysis of the previous section, but we allow /nl(^) to be an arbi- 
trary function of scale, adopting the ansatz in Eq. (2.3). We still assume homogeneity, so 
/nl(&) = /nl(^). The bispectrum is given by 

B <j ,(k 1 ,k 2 , h) = 2[h L (k 1 )P 4> (k 2 )P^(k 3 ) + perm.]. (3.15) 

Here, the triangle condition always holds, so that (for example) k\ = \k 2 + k^\. Following 
Eq. (3.7), we get 



Ab 



(k) 



D(z) ^a\M R {k) 
dfiM R {k 2 



dhkjMRik^P^h) 



/NL(A0§T?r + 2/NL(fc2) 



(3.16) 



This looks like Eq. (3.9) - but this time, /nl(^) is a function, not a constant. Thus, to 
find the derivative of Ab/b(k) with respect to the relevant parameters, we must parametrize 
/nl(&) in a way that is valid for any general form of /nl(&)- We consider the piecewise- 
constant (in wavenumber) parametrization where /nl(&) is equal to / NL in the ith wavenum- 
ber bin: 

/nl = ht(h). (3.17) 
The derivative of Ab/b(k) with respect to these / NL is: 







Off 



nl 



Ab 



(ki) 



5c 2 
D{z) 8n 2 a 2 R M R (k) 

1 



dhkfMRikJPtih) / d^M R {k 2 )P 4) {k 2 )+ (3.18) 



+2 / dhkfMRik^P^h) / dnM R {k 2 ) 

J k2£kj J 

where 5ij is the Kronecker delta function. Note that the last integral over k 2 only goes 
over the jlh wavenumber bin. 

This derivative can be rewritten more concisely as 







NL 



Ab 



25 c 
D(z) 



<5ijJi(fc) + P 2 \k) 



(3.19) 



The functions T\ and T 2 are defined as in Eqs. (3.11) and (3.12), except that the superscript 
in F% indicates that the integral over k 2 is to be executed only over the jth wavenumber 
bin. 
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4. Forecasted measurements of the scale-dependent nongaussianity 



4.1 Fisher Matrix Analysis 

With an expression for d/df^ L [(Ab/b)(ki)] in hand (Eq. (3.19)), we can calculate the Fisher 
information matrix for the parameters /i L that describe the piecewise-constant /nl(^)- 
The Fisher matrix, in turn, allows us to forecast the extent to which the scale-dependent 
non-Gaussianity could be measured in future galaxy surveys. 

We consider measurements of the power spectrum Ph(k) of dark matter halos (galaxies 
or clusters, for example) averaged over thin spherical shells in /c-space. The variance of 
Ph(k) = Ph in each shell is [60] 

_ 2 2P ft 2 /l +n p h \ 2 _ ( 2 <irP h ) 2 /1+niV 2 



h l^shell ^survey \ n Ph ) k 2 dk Vgurvey \ n Pfi ) 

where Ksheii = 47rfc 2 dk / '(27r) 3 is the volume of the shell in Fourier space (we are ignoring red- 
shift distortion effects for simplicity here). Therefore, the Fisher matrix for measurements 
of P h (k,z) is [61] 



~vi " fan 



dP h (k,z m ) dP h {k,z m ) 1 k 2 dk 



dpi dpj 



Ph(k,z m ) + - 
n 



2 (2tt) 2 



(4.2) 



where V m is the comoving volume of the m-th redshift bin, each redshift bin is centered 
on z m , and we have summed over all redshift bins. We adopt k m \ n = 10~ 4 /i" 1 Mpc, 
and we choose k max as a function of z so that (j(Tr/(2k max ),z) = 0.5 [62], which leads to 
fc max (z = 0) ~ 0.1/iMpc -1 . Finally, pi are the parameters of interest; in our case, these 
are the f^ L . 

We assume a flat universe and a fiducial model of zero non-Gaussianity: /nl(&) = 
= /jvjl- We include six cosmological parameters in our Fisher matrix aside from the 
/^ L : Hubble's constant Hq; physical dark matter and baryon densities Q c dmh 2 and Q^h 2 ; 
equation of state of dark energy w; the log of the scalar amplitude of the matter power 
spectrum, log A s ; and the spectral index of the matter power spectrum, n s . Fiducial values 
of these parameters correspond to their best-fit WMAP7 values [63]. We also added the 
forecasted cosmological parameter constraints from the CMB experiment Planck by adding 
its Fisher matrix as a prior (W. Hu, private communication). Note that the CMB prior 
does not include CMB constraints on non-Gaussianity; the CMB constraints on /nl(^) will 
be separately studied in a future work. Finally, in addition to the cosmological parameters 
and the f^ L , we include five Gaussian bias parameters in our Fisher matrix - one bo(z) 
for each redshift bin. The fiducial values of these parameters are set by the relations 
b ( z = 0) = 2.2, and b (z) = b (z = 0)/D(z). 

We already have the derivatives of b(k) with respect to each of the /^ L , so the derivative 
of Ph(k) with respect to the is just 

?lf>^m b{k)p ^ (k); (4 . 3) 
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-fmat(^) is the ACDM matter power spectrum, easily obtained from a numerical code such 
as CAMB. Since we only consider information from large scales (k < k max ~ 0.1 /iMpc -1 ), 
we do not model the small amount of nonlinearity present at the high-/c end of these scales. 

We assume a future survey covering one-quarter of the sky (about 10,000 square de- 
grees) out to z = 1, and find constraints for a set of 20 f^ L uniformly spaced in log A; in 
the range 10 -4 < fc/(/iMpc _1 ) < 1, with a smoothing scale of M smoot h = 1O 14 M . Fig. 1 
shows the resulting unmarginalized (left panel) and marginalized (right panel) constraints 
on the parameters f^ L . For both sets of constraints, we first marginalized over the other 
cosmological parameters. 3 The have most of their degeneracy among themselves; a 
plot showing the fully unmarginalized constraints on the f^ L would not look much dif- 
ferent than the left panel of Fig. 1. Note that, while some of the have support at 
k > k maiX (z = 1) « 0.2frMpc _1 , we only use information about those (and other) pa- 
rameters coming from k < /c max - The constraints vary considerably as a function of the 
k at which these parameters are defined. The best-constrained corresponds to the 
1Q-0.8 ^ io~ ' 6 bin, and it has an estimated unmarginalized error of o"(/nl) = 7-3; for 
comparison, the worst-constrained f^ L , which corresponds to the largest scale (smallest k) 
bin, has an unmarginalized error well over one billion. 

As expected, the marginalized constraints for the best-constrained parameters are 
much weaker than the unmarginalized constraints - even the best-measured has an 
estimated marginalized error of 6 x 10 2 . In general, dependence of the constraints on the 
value of k is determined by two competing factors: as k increases, there is a larger number 
of modes, each with a smaller signal (given by the smaller nongaussian bias A6). The 
best-constrained k is also affected by the fact that only information out to k = k max = 
0.1/iMpc" 1 is assumed from the galaxy survey. In particular, we have checked that if we 
unrealistically assume information to be available at all k (instead of at k < fc max ) without 
modeling the nonlinearities, the unmarginalized constraints on f^ L improve monotonically 
with increasing k. Therefore, the raw signal-to-noise ratio in f^ L increases with k. To 
further demonstrate the effect of the choice of k max (z), we also plotted the errors obtained 
with the condition o"(7r/(2/c max ), z) = 0.15, which yields fc max (z = 0) « 0.03. 

The smoothing mass scale chosen for this analysis (see Eq. (2.5)) has a small but 
noticeable effect on the constraints yielded. Figure 2 shows that, in the case of the un- 
marginalized errors, the k at which non-Gaussianity is best constrained decreases as the 
smoothing mass scale increases. (The behavior of the marginalized errors is more com- 
plicated due to correlations in errors between neighboring /n L -) Since the mass scale is 
proportional to the physical scale (to the third power), this means that best-constrained 
k decreases with increasing smoothing scale R, which is exactly what we should expect. 
We remind the reader that while a survey filtered at some scale M smoot h contains ob- 
jects roughly more massive than this scale, in practice the near-exponentially falling mass 
function implies that the number density is dominated with M ~ M smoo th halos. 

3 Using six cosmological parameters along with five bo(z) and 20 /Jj L led us into some issues with floating- 
point errors and numerical precision. The 31 x 31 Fisher matrix we obtained was rather ill-conditioned and 
difficult to invert reliably using 64-bit precision; we were eventually forced to move to 128-bit precision in 
order to accurately marginalize over the cosmological parameters. 
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Figure 1: Estimated unmarginalized (left panel) and marginalized (right panel) constraints on 
piecewise-constant parameters /^j L assuming a future galaxy survey covering one-quarter of the sky 
out to z = 1, with average number density of 2x 10 -4 gal/Mpc 3 . For comparison, the green line is the 
constraint found for a constant /nl using the same survey assumptions, and the red histograms are 
the constraints found with a lower /c max (see text for details) . While the individual parameters 
are poorly constrained as expected, their few best linear combinations - the principal components 
- are well measured; see the next section and text for details. 
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Figure 2: Estimated constraints obtained from future surveys with the same parameters as the 
previous figure at different mass smoothing scales M smoo th (labeled as M in the legend). In other 
words, these are errors for a survey with halos of M > M smoo th- 



5. Projection and Principal Components 
5.1 Constraining other /nl(&) models 

Once the Fisher matrix F has been obtained for the set of parameters /L, it is quite simple 
to find the best possible constraints on the that could be obtained from a future galaxy 
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redshift survey. By projecting this Fisher matrix into another basis (see Appendix A), it 
is also possible to find the constraints on any arbitrary /nl(&) without calculating a new 
Fisher matrix from scratch. A trivial example can be found in Appendix A, where we find 
that the estimated error on a constant /nlj assuming the same future survey as in the 
previous section, is c(/nl) — 2.1. (Note that this forecasted constraint is on a par with 
the error expected from Planck, where <t(/nl) ~ 5.) 

For another, scale-dependent example, consider the simple form of non-Gaussianity 
analogous to the conventional parameterization of the power spectrum 

/ k \ " NG 

/nl(*0 = /nl ( Y ) > ( 5J ) 

where fc* is an arbitrary fixed parameter, leaving and h-ng as the parameters of interest 
in this model, (fc* is generally chosen to minimize degeneracy between f^ L and hng for 
the observable of interest. We have set k* = 0.165ft, Mpc -1 , close to the optimal value in 
our case; in CMB analysis, the optimal value is lower, around O.O6/1 Mpc -1 .) The partial 
derivatives of our basis of / NL with respect to these parameters are: 

^% = (-V NG ; (5.2) 

0^ = /nl UJ s UJ • (5 - 3) 

Starting in a basis of 20 / NL evenly spaced in log k, we project down to a basis of / NL and 
ting m order to forecast constraints on the two new parameters from a survey covering one- 
quarter of the sky out to z = 1. We are using the same limits of integration as in Section 
4.1, along with the fiducial values / NL = 50 and 7?ng — 0. The forecasted constraints on 
these parameters, marginalized over each other, are cr/* = 1.7 and o" n . NG = 0-58. Despite 
a superficial similarity between this model and the model used by Sefusatti et al. in [23], 
the two models are quite different, and our results cannot be compared. The model used 
in [23] is a function of three arguments, k\,k2, and ^3: 

f K\ nNG 

fNL(h,k 2 ,k 3 ) = f^ L ^-j , (5.4) 

where K = (fci^^s) 1 ^ 3 . This leads to a bispectrum of the form found in Eq. (3.8), but 
with /nl(^1) &2) ^3) in place of /nli whereas our bispectrum is of the less-factorizable form 
Eq. (2.7). 

Another example we consider is the form of non-Gaussianity in which the running on 
/nl itself has running; that is, the case in which ting is a function of k. A simple case of 
this would be /nl of the form 4 

hL(k) = e AkB . (5.5) 

Projecting the Fisher matrix down from the original basis /^ L to the parameters A and B, 
with fiducial values of A = log 50 and B = 0, we obtain forecasted constraints of a a = 1.0 
and as = 0.15. (In this case, the survey characteristics and bounds of integration are the 
same as in the previous example.) 

4 Analogous parameterization for the power spectrum and its motivations are discussed in [64]. 
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Figure 3: The first four principal components of /nlCO- The PCs, e^\k), are eigenvectors of 
the Fisher matrix for the /n L , and are ordered from the best-measured one (j = 0) to the worst- 
measured one (j = 19) for the assumed fiducial survey. 

5.2 Principal components and relation to local and equilateral models 

We now represent a general function /nl(^) i n terms of principal components (PCs). In this 
approach, the data determine which particular modes of /nl(&) are best ° r worst measured. 
The PCs also constitute a useful form of data compression, so that one can keep only a 
few of the best-measured modes to make inferences about the function /nl(&)- Finally, the 
PCs will also enable us to measure the degree of similarity between our scale-dependent 
ansatz and the local and equilateral forms of non-Gaussianity. 

It is rather straightforward to start from the covariance matrix for the piecewise con- 
stant parameters fist, an d obtain the PCs of /nl(&)- The PCs are weights in wavenumber 
with amplitudes that are uncorrelated by construction, and they are ordered from the 
best-measured (i = 0) to the worst-measured (i = 19) for the assumed fiducial survey. 
The construction of the PCs is described in Appendix B. A few of these PCs of /nl(&) 
are shown in Fig. 3. For example, the best-measured PC has most of its weight around 
k = 10- - 4 /iMpc _1 , which agrees with sensitivities of piecewise-constant parameters shown 
in Fig. 1. The sensitivity is not greatest at the largest value of k (1 /iMpc -1 ) because we 
assumed cosmological information from k < k mSLX = 0.1/iMpc -1 . We checked that in- 
formation available at a higher A; max would shift the "sweet spot" of sensitivity to higher 
wavenumbers. 
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Figure 4: RMS error on each principal component, along with the cumulative error. 



The error in the best-measured PC is 4.8; however, the error in the next-best measured 
PCs are 18.3 and 27.4, and the accuracy rapidly drops off from there. Thus, the first three 
or four PCs should be enough for any conceivable application. The error in each PC is 
plotted on a logarithmic scale in figure 4, along with the cumulative error a cum , which is 
defined as 

(5.6) 



Each PC e^(Jfe) has its own associated bispectrum (see Eq. (2.7)): 

B^(k 1 ,k 2 ,h) = 2[e^(k 1 )P(k 2 )P(k 3 ) + e^(k 3 )P(k 1 )P(k 2 ) + e^(k 2 )P(k 3 )P(k 1 )]. (5.7) 

(As always, k\,k 2 , and k 3 have a triangle relation: k 3 = \k 2 — k\\.) We would like to test the 
similarity of these bispectra to those that have already been discussed in the literature. We 
can do this by using a distance measure between bispectra, defined by 'cosines' developed 
in [40]. A cosine near unity implies that the two bispectra have very similar shapes, and a 
cosine near zero implies the opposite. The cosine is defined as 

P> B 

cos(B 1 ,B 2 )= , 1-2 (5.8) 
y /(B 1 -B 1 )(B 2 -B 2 y 

where the inner product between two bispectra, B\ ■ B 2 , is [23] 

D D B 1 (k 1 ,k 2 ,k 3 )B 2 (k 1 ,k 2 ,k 3 ) 

The (Gaussian) variance of the bispectrum is 

A 2 B(h, k 2 , k 3 ) = -L P {k l )P{k 2 )P{k^) ~ j^(hk 2 k 3 y 3 , (5.10) 
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Local cosine 


Equilateral cosine 




0.669 


0.074 




0.040 


0.000 




0.099 


0.030 


B (3) 


0.189 


0.037 



Table 1: Cosines of the first four principal-component derived bispectra with the local bispectrum 
and the equilateral bispectrum. A cosine near unity implies that the two bispectra have very similar 
shapes, and a cosine near zero implies the opposite. Note that the zcroth PC, which is by far the 
best measured (see Fig. 4) , has a much larger overlap with the local model than with the equilateral, 
as expected. 



where Nt is the number of distinct triangular configurations of ^2,3, and P(k) ~ k~ 3 is 
the primordial curvature perturbation power spectrum. (The overall constant is irrelevant, 
since it cancels out in Eq. (5.8).) 

We first compare our bispectra Eq. (5.7) to the local model with a constant /nl> whose 
bispectrum is (see Eqs. (2.1) and (3.8)) 

111 

Bi oca x{k 1 ,k 2 ,k 3 ) oc pp+pp+pp- (5.11) 

12 13 2 3 

Most of the power of -Biocal is in so-called "squeezed" triangles, in which one side is much 
smaller than the other two (comparable) sides, k\ « k 2 ~ ^3. 

Another form for the bispectrum much discussed in the literature is the "equilateral" 
bispectrum 



1 

(feAqfo) 2 locau Li z ' aj ' kxk^k^ ' k^kfk'^ 



B eqni (ki,k 2 ,h) = - n _ - - ^ -Bi oca x{ki,k2,kz)+ T , 2 ,^ + , , 2 , .^ +permutations. (5.12) 



In contrast with B\ ocs \, most of the power of -B equ i is in triangles where k\ ~ k 2 ~ ^3; hence 
the name "equilateral". 

Table 1 lists the cosines of the first few principal-component derived bispectra with 
the local bispectrum and the equilateral bispectrum. The form of Eq. (5.7) suggests that 
the PC-derived bispectra B^J will have more in common with the local bispectrum than 
the equilateral one. However, it is initially conceivable that some e^\k) might exist which 
would yield a bispectrum of the form in Eq. (5.12) when substituted into Eq. (5.7) - but 
in Appendix C, we prove that no such function exists. Thus, the only guarantees for the 
cosines of the are that the cosine of B^ - the bispectrum corresponding to the best- 
measured PC - will be large with the local model, and that none of the B^ have a very 
large cosine with the equilateral model. We expect the former because our model looks 
like the local model; we expect the latter because of the proof in Appendix C. Table 1 
bears out this expectation. The small cosines with the equilateral form of non-Gaussianity 
are also unsurprising because equilateral non-Gaussianity is expected to have a strongly 
suppressed signal in the non-Gaussian halo bias [57]. 
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6. Conclusions 



In this paper we have suggested a new phenomenological model of primordial nongaus- 
sianity by generalizing the local model (parametrized with a constant parameter /nl) to 
a scale-dependent, non-local class of models. There are multiple ways to do this, and our 
choice was to write the Newtonian potential as 

$(x) = <j> G {x) + /nl(x) * {^{xf - (4>g(x) 2 }), (6.1) 

where the convolution in real space corresponds to multiplication in fc-space, featuring 
an arbitrary function /nl(&). Explicit calculations show that such a form of the scale 
dependent /nl is borne out in inflationary models [11, 14, 37, 38, 39]. 

We calculated the bispectrum and bias of dark matter halos in this class of mod- 
els, following the formalism valid for high peaks [58, 59]. We then specialized in the 
piecewise-constant (in wavenumber) parametrization of /nl(&) which, for the case of nar- 
row enough £>bins, recovers any arbitrary function. We used forecasted constraints from 
an intermediate-future galaxy survey to calculate errors on individual parameters / NL (see 
Fig. 1) and briefly studied dependence on the smoothing scale (Fig. 2). 

We further calculated the principal components of /nl(&)> an d thus identified the 
best-measured configurations (in wavenumber) of this function (see Fig. 3). While the 
sensitivity increases with increasing k, restricting the survey information to scales where 
linear perturbation theory is valid imposes a "sweet spot" in sensitivity of k ~ 0.1/tMpc - . 
We then calculated the overlap of the best-measured principal components with two fa- 
miliar classes of non-Gaussian models: local (/nl = const) and equilateral models, using 
a cosine measure between the bispectra suggested in [40]. We found the expected result: 
the best measured component overlaps much more with the local model (which our model 
generalizes) than with the equilateral one. 

One immediate utility of our results is an easy adaptation to specific models of non- 
Gaussianity predicted by classes of inflationary models. If one wants to forecast the ac- 
curacy with which parameters of a specific model of /NL(&)-style non-Gaussianity will be 
measured, neither the halo bias nor the Fisher matrix needs to be calculated from scratch. 
Instead, our formalism makes it possible to obtain these forecasts by performing a simple 
linear projection to our piecewise-constant model; this procedure is described in Appendix 
A and illustrated with a few examples. 

In future investigations, it will be interesting to consider specific inflationary models, 
projecting down to specific forms for /nl(&)- It will also be important to test how well the 
observable effects of scale-dependent non-Gaussianity, studied here using the theoretical 
ansatz from Eq. (3.2), agree with numerical simulations; the first such investigations, for 
select specific forms of /nl(&)> are now being done [41]. Finally, it will be interesting to see 
how one can optimally select objects in the universe (i.e. their mass) to probe information 
about scale-dependence of non-Gaussianity. While in Fig. 2 we showed scaling of the best- 
determined scale of /nl(^) with the smoothing mass scale applied to the density field, a 
more complete analysis might use the Halo Occupation Distribution (HOD) approach to 
relate the content of dark matter halos to their mass. 
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A. Calculating the error on an arbitrary parametrized /nl(&) 

Projecting the constraints from an old set of parameters f^ L = /nl(^i) i} = 1> 2, . . . , iV) 
to new parameters (which we can call q; j = 1,2, ... ,M for some M) is in principle 
straightforward. The Fisher matrix in the new parameters, F new , is given by 

k,l=l 

so that 

F new = V T FV, (A.2) 

where V%j = dp l /dq 3 is the derivative matrix of old parameters with respect to new. 
Let us look at a couple of examples. Projecting to the case 

hh(k) = /nl = const (A.3) 

is particularly easy, since V is the column vector with V%\ = ^/nl/^/nl = 1- Then F^ cw is 
a 1 x 1 matrix that quantifies information on /nl> given by 



k,i 



The error on /nl is of course given simply by <t(/nl) = l/\/^ii CW 



Another example is given by the function 



k \ nNG 

hL(k)=[ Tr ) , (A.5) 



with two parameters, fco and n^G- Then one can show that (labeling k$ = q\ and ?t,ng = Q2)' 

Va = - *T UJ ' (A - 6) 

{£T- (A - 7) 

Then, using Eq. (A.2), one can simply obtain the 2x2 Fisher matrix in ko and ting- 
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B. Principal Components of /nl(&) 

We now show how to decompose the measurement of /nl(&) m principal components, 
which are essentially the eigenmodes of the covariance matrix for the aforementioned pa- 
rameters /nl(&i). This method has been widely used in cosmology, including applications 
to parametrizing and describing dark energy [65, 66]. It allows us to order the best-to- worst 
measured weights in wavenumber of the function /nl(&)- 

Let the function /nl(&) be described in terms of piecewise constant parameters ffc L = 
hhih), where 

N 

/ NL (fe) = 5>ei(fc). (b.i) 

i=i 

Here, @(k) = [H(k — k\ owcr ) — H(k — A;" ppcr )] is the top-hat function of unit height over 
the ith wavenumber bin, and we assume a total of N bins. k\ owcv and fc" pper are the 
wavenumber bin boundaries, and H is the Heaviside step function. We have effectively 
expanded the function around the zero value, though this is not crucial: the left-hand side 
could be /nl(&) — /nl(^)> f° r an y fiducial f^{k), and the formalism still follows. 

The Fisher matrix F is the inverse covariance matrix in the original piecewise-constant 
parameters pt, so that F^ 1 = (piPj) — (pi)(pj)- We first diagonalize the Fisher matrix F: 

F = W T DW, (B.2) 

where D is diagonal and W is some orthogonal matrix. The vector of uncorrelated param- 
eters, q, is related to the vector of original parameters p via 

q = Wp, (B.3) 

and it is easy to check that the q are uncorrelated; that is, (qq r ) = D^ 1 . The rows of W 
are therefore the new parameters. 

Thus, to calculate the principal components: 

1. Obtain the full Fisher matrix for iV parameters pi, plus the cosmological parameters 
Q b h 2 ,Q C DMh 2 ,H ,w, log A s , and n s . 

2. Marginalize over the cosmological parameters by inverting this larger Fisher matrix, 
taking the N x N submatrix, then inverting back to get the Fisher matrix of the pi\ 
we call this Fisher matrix F 

3. Diagonalize F as in Eq. (B.2) 

4. The rows of W are the principal components. More precisely, q a = ^ W a iPi, and q a 
are the PCs. 

Let us now change notation slightly (to agree with the commonly used one, e.g. [65]), 
and define the shape of the a-th principal component in i-th. redshift bin as oq \ so that 
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a ( a ) = W a i. Then we can represent the a-th principal component, e^ a \k), in terms of the 
original parameters pi as 5 

N 

e^(k) = Y,^ ) p i Q i (k). (B.4) 

i=l 

The PCs are obviously uncorrected, and their eigenvalues A a , so that 

(eVeV)=j24 a) <*?<PiPj) = ¥- (B-5) 

where, recall, A a = D aa . 

Finally, let us calculate the coefficients in the expansion in principal components 
of an arbitrary /nl(&) 

N 
a=l 

Let coefficients /^l describe /nl(&) in our original basis, so that /nl(&) = const = 
Si fhhPi®i(k)-> /nl being left arbitrary for now. Then, taking the expectation value 
of the product with e^ b \ we get 



N \ / AT N 



=1 



AT 



so that 



N 

„(<*) C 



For example, in the simplest case of constant /nl(&), where / NL = const = /nl, the 
coefficients of the principal components in the expansion of /nl(&) are 



Ca = K /NL E ^ ( fOT = ^ NL = COnSt )' ( R10 ) 



C. Generalized local ansatz does not recover the equilateral case 

Here, we prove that our ansatz cannot perfectly mimic the equilateral bispectrum for any 
choice of /nl(&)- The generalized local form of the bispectrum that we considered in this 
paper is 



f ) 

Bgener(fci, k 2 , fo) = 2[fah(ki) P (k 2 ) P (k 3 ) + permutations] oc ^ 3 3 * + perm. (C.l) 

fc 2 fc 3 



3 This is basically the continuous version of the relation q a — J2i WaiPi- 
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The equilateral bispectrum is 



B equi (ki, k 2 , k 3 ) oc 



k\ krj k^ 



+ perm. 



(k 2 k 1 k 3 y 



h 3 h 3 



+ perm. 



(C.2) 



The claim is that there is no /nl(^) such that -B gen er = -Bequi for all k\,k 2 ,k 3 - To show 
this, we define a new function h(k) = fNL(k) + 1. If there is some /nl(^) such that 



^Bgcner = -Bequi, then we have: 



hjh) 
h 3 h 3 



+ perm, oc 



k\ k<2 k^ 



+ perm. 



(k 2 kik 3 )< 



We can go from a proportionality to an equality by defining a new function g{k) that 
is simply h(k) with the appropriate constant out in front. Next, multiply both sides by 
kfk^k^ to get 



klgikx) + klg{k 2 ) + klg(k 3 ) = \k x k\ + k 2 k\ + perm.] - 2k x k 2 k z . 



(C.3) 



Each term on the left-hand side is dependent on only one of k\,k 2 , or ^3. However, every 
term on the right-hand side depends on at least two different k; thus, there is no g{k) that 
can satisfy this relation. 

Alternatively, consider the case where k\ = k 2 = k% = k. Then (C.3) becomes 



which means that 



3k 3 g(k) = Ak 3 



g(k) = 4/3. 



This answer is wholly independent of k, so this value of g{k) must be true for all k. But 
this solution for g{k) is clearly incorrect in the general case where k\ 7^ k 2 7^ therefore, 
no such g(k) can exist. 

While this proves that there is no /niX^) that yields an exact equality between our 
ansatz and the equilateral bispectrum, the question of an approximate equality remains. 
Such solutions for /nl(^) certainly exist for narrow ranges of k. For example, /nl(^) = 
5(k—k*), where 5(k) is the Dirac delta function, yields a bispectrum that is larger for exactly 
one equilateral triangle - the triangle where ^2,3 = k* - than it is for any squeezed triangle. 
However, no /nl(^) exists that yields a bispectrum which favors equilateral triangles over 
squeezed triangles for all k. It is straightforward but tedious to prove this fact, and the 
details of the proof are beyond the scope of this paper. 
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